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ABSTRACT 

Context. Solar tornados are dynamical, conspicuously helical magnetic structures mainly observed as a prominence 
activity. 

Aims. We investigate and propose a triggering mechanism for the solar tornado observed in a prominence cavity by 
SDO/AIA on September 25, 2011. 

Methods. High-cadence EUV images from the SDO/AIA and the Ahead spacecraft of STEREO/EUVI are used to 
correlate three flares in the neighbouring active-region (NOAA 11303), and their EUV waves, with the dynamical 
developments of the tornado. The timings of the flares and EUV waves observed on-disk in 195A are analyzed in 
relation to the tornado activities observed at the limb in 171 A. 

Results. Each of the three flares and its related EUV wave occurred within 10 hours of the onset of the tornado. They 
have an observed causal relationship with the commencement of activity in the prominence where the tornado develops. 
Tornado-like rotations along the side of the prominence start after the second flare. The prominence cavity expands 
with acceleration of tornado motion after the third flare. 

Conclusions. Flares in the neighbouring active region may have affected the cavity prominence system and triggered the 
solar tornado. A plausible mechanism is that the active-region coronal field contracted by the 'Hudson effect' due to the 
loss of magnetic energy as flares. Subsequently the cavity expanded by its magnetic pressure to fill the surrounding low 
corona. We suggest that the tornado is the dynamical response of the helical prominence field to the cavity expansion. 
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1. Introduction 

Prominences consist of relatively large, cool and over 
dense plasma seen in the lower corona above the solar 
limb (|Martinlll998l : lTandberg-Hanssedll995t iMackav et all 
l2010l ). Their structure and composition are exceedingly 
complicated. The pla sma mainly resides in highly tan - 
gled magnetic fields ([van Ballegooijen fc Cranmerl 120101 ). 
Coronal cavities are often observ ed to have cooler 
prominence plasm a at their bases (Hudson et al.l [l999; 
iGibson et alTl2006l : iRegnier et all |20n|). The cavity is a 
region of relat i vely l o w density, high temp er ature plasma 



egion 01 r elatively lo w density, nign temperatu re plasma 
Gibson et all 120101 lHabbal et all |2010[ ) . iBerger et al.l 
( 201l[ ) have proposed that the prominence and its cavity 
are a form of magneto- thermal convective structure, macro- 
scopically stabl e but internally in a constant state of ubiq- 
uitous motions (|Low et alJl2012allbh . 

Viewed at the limb, quiescent prominences mostly 
appear as curtain s of vertical thread-like structures 
(|Berger et al.| [2008). Occasionally they look like to rnados 
with r otations along, or of, t heir magnetic structures ([Pettitl 
119431 : iLiggett fc Zirinl 119841 ). At first glance it is not ob- 
vious why tornados do not erupt but continue rotating 
for several hours before quietening down. Recently, im- 
ages from Atmospheric Imaging Assembly (AIA) on the 
Solar Dynamic Observatory (SDO), have allowed the study 
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of sola r tornados (|Li et alJ[20ll ISu et al.|[20ll iLiu et al.l 
l2012af ). The driving mechanism may be a coupling and 
exp ansion of a tw isted flux rope into the coronal cav- 
ity ([Li et al.l l2012f ) and/or related to photospheric vor- 
tices at the footpoint of the tornado (I Attie et all [2009; 
IWedemever-Bohm et al.ll20H ISu et al]l2Q12~ ~~ 

This letter is a study of the impressi ve solar t o rnado 
observed by AIA on 25 September 2011. iLi et al.l (|2012[ ) 
described the formation and disappearance of the related 
prominence by analysing SDO data from 24-26 September. 
The evolution of the tornado was attributed to the expan- 
sion of helical structures into the cavity. But the reason 
behind the expansion remained an open question. 

Here we analyze observations on the mechanism leading 
to the expansion of the prominence. We show using Solar 
TErrestrial RElations Observatory (STEREO) Extreme 
UltraViolet Imager (EUVI) observations that three strong 
flares in a neighbouring active region coincided with phases 
of the tornado activation. Each flare was associated with a 
Coronal Mass Ejection (CME) and EUV wave. Particularly, 
after the third flare, there was slow but significant expan- 
sion of the overlying prominence cavity together with more 
rapid rotation at the top of the prominence. We describe the 
observations in Section 2. The evolution of the tornado, its 
relationship to the solar flares and EUV waves are described 
in Section 3. In Section 4, we summarize our observations 
and speculate on the link between the flares and the cavity 
expansion. 
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Fig. 1. The region 
around a tornado on 
25 September 2011: 
a) EUVI-A 195A; b) 
AIA 17lA images. 
On the EUVI image, 
the black dashed line 
shows the position of 
the AIA limb, and the 
long diagonal arrow is 
the epipolar line for 
the prominence posi- 
tion marked with '+' 
on the AIA image. 
The orange dashed 
line outlines the po- 
sition of the south- 
ern edge of the AR 
corona. 



2. Observations 

The solar tornado on 25 September 2011 was observed on 
the south west limb by SDO/AIA (jLemen et al.ll2012l ). In 
images from EUVI on the Ahead spacecraft of STEREO 
([Howard et al.ll2008l ), it appeared on the solar disk around 
15° E. The separation angle between SDO and STEREO-A 
was 103°. We study the tornado by combining observations 
from the two directions. 

AIA takes high spatial resolution (0.6" pixel -1 ) full disk 
images with a cadence of 12 s. For the analysis, we select 
images from the 171 A channel which is centered on the Fe IX 
line formed around 0.63 MK. Images from this channel show 
both hot, bright emitting and cold, dense absorbing parts 
of the prominence, as well as the dark, low density cavity 
and the surrounding coronal loops. To enhance the visibil- 
ity of the cavity and faint loops, we removed an average 
background from the 171 A images by taking the median of 
two months data - September and October, 2011. 

During the period studied, 01:30-13:30 UT, SDO was 
in eclipse from 06:02-07:13 UT. This gap is covered b y 
the SWAP instrument on PROBA-2 ([Halain et al.ll2010l ). 
SWAP provides a full Sun 174 A image every 2-3 min with 
a spatial resolution of 3. 16" pixel -1 . The images are not as 
detailed as AIA images but are essential for checking the 
behaviour of the prominence during the data gap. 

The EUVI-A 195A images have a time cadence of 
5 min and resolution of 1.6" pixel -1 . The 195 A emission 
is mainly Fexn formed at ~1.2 MK. To confirm the promi- 
nence position in the EUVI images, we obtained the three- 
dimensional coordinate s of the prominen ce with the rou- 
tine SCCLMEASURE ([Thompson! I2009D available in the 
SolarSoft library. 

2.1. Overview 

In Fig. [1] we show the EUVI-A 195A and AIA 171A images 
of the filament /prominence and the active region NOAA 
11303 (marked with an arrow) as seen from the two angles. 
The position of the cross on the stem of the prominence 
lies along the epipolar line drawn as a long black arrow 
on the EUVI image. The dark lane crossed by the epipolar 
line is the filament channel and the tornado is the bright 
region that coincides with the epipolar line and the fila- 
ment channel, indicated with a black arrow. The core of 



the active region is separated by 300" from the filament 
channel. However, the southern edge of the active region 
corona (outlined by orange dashed line) reaches to within 
50" of the filament channel. 

2.2. The three flares 

There were flares at approximately 02:45, 07:00, 09:40 UT 
from the nearby active region. They were all associated 
with CMEs and EUV waves. The third was GOES class 
Ml. 4 according to the NOAA flare catalogues The GOES 
class of the first two are not given in the catalogue. We 
therefore checked the hard X-ray quicklook images from the 
Reuven Ramaty High Energy Solar Spectroscopic Imager 
(RHESSI). Unfortunately they both had their peaks during 
RHESSI data gaps. The first RHESSI image in the flares' 
decay phase showed that these flares were the brightest 
RHESSI sources. At the time of the first and second flares 
the GOES 1-8 A fluxes reached the M4.4 and M1.0 level 
respectively, so it is possible that these too were M-class 
flares. 

The flares' onsets and related EUV waves can be seen 
in the 195A running ratio images shown in Fig. [2j The run- 
ning ratio images are the log of the intensity at the time 
shown on the images divided by the intensity 5 min earlier. 
Because this reflects relative changes the faint EUV waves 
show up against bright inactive active regions. The faint 
EUV wavefronts sweep over the filament channel about half 
an hour after the flare (bottom row). The small white ar- 
rows point to the tornado site where there was increased 
activity along the filament channel after the flares. 

When we take a time series along the orange dashed 
line drawn on Fig. [2^, we are able to see the relation be- 
tween the EUV waves and the filament activity (Fig. [3]). In 
this time series image, the waves associated with the three 
flares are marked by Fl, F2 and F3 and a red arrow points 
to the filament activity. We note that the filament activity 
associated with F2 increased before the arrival of the EUV 
wave. We also notice that each EUV wave front was fol- 
lowed by a series of oscillations with a period of 20 min at 
the edge of the active region and along the filament chan- 
nel. Similar post-EUV cavity oscillations, seen at the limb 
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Fig. 2. EUVI-A 195A running ratio images showing the onset of the three flares (a,b,c) and their associated EUV waves 
(d,e,f). In (a) the epipolar line shown in Fig. [TJi is drawn as a long white arrow, and the solar limb is drawn as a dashed 
white line. The small white arrow in all images points to the tornado site. The dot-dashed/dashed lines in (b) mark the 
positions of time series shown in Fig. |4j and the orange dashed line along the EUV wave front in (e) the position of the 
time series in Fig. [3j 
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Fig. 3. EUVI 195A running ratio space time image along 
the dashed line in Fig. [2^. Three flares are labelled as Fl, 
F2 and F3. The red arrow points to the filament. 



in AIA images, have been reported by iLiu et al.l (|2012bl ). 
EUVI images have a much lower cadence than AIA, and 
may not be resolving the wave trains. 

As mentioned above, the EUV wave reached the fila- 
ment channel after the filament activation. To investigate 
whether the activation could have been triggered by the 
flare, we took time series through the flare sites (vertical line 
Fig. [2}d) and along the filament channel across the activated 
filament (horizontal line in Fig. [2Jd). We observe that there 
are perturbations in the filament plasma immediately after 
each flare (yellow arrows in Fig. HJd). Large-scale changes 
along the filament channel started after F2 and increased 
further after the F3. 

The reaction of the prominence to the flares is illus- 
trated in Fig. The positions of the three time series are 
drawn on Fig. [6j an AIA image of the active region and 
tornado. The oscillations of the prominence at the time of 
Fl (Fig. [T^i) are consistent with a flare trigger from the 



north since it moved first towards the pole. The wobbling 
in the prominence after the flare is also clearly visible in 
the movie6a (Fig. [6]). 

The second flare, F2, was at -06:45 UT when SDO was 
in eclipse. Since this phase is important for understand- 
ing the tornado, we checked the SWAP data to see if the 
prominence activation started before F2. The SWAP data 
are shown in the movie attached with Fig. [5j Although the 
SWAP images are not as sharp as AIA images, they show 
quite clearly that increased prominence activity occurred 
after F2 and that the first AIA images after eclipse caught 
the expansion of the flare loops and the rapid growth in 
prominence activity. The middle image in Fig. [TJd, taken 
along the dashed line in Fig. [6j shows an arm of promi- 
nence plasma reaching out in the direction of the active 
region. There is also a very slight expansion of the cavity 
boundary towards the active region. 

The biggest change in the cavity occurred after F3 when 
the bright and complex top part of the prominence started 
to rotate faster. The movie6b attached to Fig. [6j shows the 
evolution of the tornado and flare from F2 to a couple of 
hours after F3. The time series in Fig. taken along the 
top white arrow in Fig.[6j shows the main features. The flare 
erupted at 09:20 UT. The active region loops started ex- 
panding at 09:40:48 UT. At 09:50:13 UT there was a sharp 
contraction, approximately 10 Mm, of the observed cavity. 
The cavity started expanding back towards the active re- 
gion after —10:00 UT. It grows from 140 Mm to 167 Mm in 
3.5 hours. At the same time, the prominence grew in height 
and the bright and complex head of the prominence started 
to rotate faster, forming the tornado. The rotations at the 
top of the prominence/tornado last ed —3.5 hours with a 
speed from 55 km s" 1 to 95 km s" 1 ([Li et al.ll2012l ). 
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Fig. 4. EUVI 195A running difference time series along (a) 
the vertical and (b) the horizontal line in Fig. [2}d. The red 
dashed lines are drawn at the flare Fl, F2, and F3 times. 
In b) the yellow arrows point to the changes in filament 
structure after the flares. 
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Fig. 5. SWAP 174 A intensity image from 25 September 
2011. The arrows show the positions of the cavity and the 
active region (AR). This frame is taken from the movie- 
"MOVIE5.mp4". 



3. Discussion 

A solar tornado was observed by AIA and EU VI-A on 
25 Sep tember 2011. The tornado was reported bv lLi et all 
(2012), who attributed it to the growth of a helical promi- 
nence system. The reason for the prominence growth re- 
mained unanswered. After careful observation of the event 
using AIA full disk images, we noticed that AR 11303, to 
the north of the prominence cavity was flaring and sus- 
pected it of having influenced the prominence in some way. 
We observed three flares, and they were all associated with 
EUV waves that swept over the prominence cavity. Each 
time the cavity was buffeted by the wave, the prominence 
plasma became more active. After the third flare a tornado 
had developed at the top of the prominence. 

After the first flare, there were clear oscillations along 
the prominence stem. The second flare caused plasma to 
move along arm-like extensions projecting out from the 
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Fig. 6. AIA 171 A intensity image. The diagonal long ar- 
rows show the positions of the time series images through 
the prominence/tornado, cavity, active region hot loop (AR 
loop) and active region (AR) shown in Fig. [71 This is the 
frame from the movie-" MOVIE6b.mp4" . "MOVIE6a.mp4" 
shows the evolution of this region at earlier times. 
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Fig. 7. AIA 171 A intensity time series along the diagonal 
arrows in Fig. [6l a) oscillations of the prominence stem 
triggered by Fl taken along bottom white arrow; b) activity 
triggered by F2 and F3 taken along the middle white arrow; 
c) tornado activity and cavity changes triggered by F2 and 
F3 taken along top arrow. The arrow at the time of the 
flare (~09.40 UT) points to the movement of the cavity 
boundary. 



stem in the direction of the active region. These extensions 
appear to subsequently rotate back towards the prominence 
before the third flare. After the third flare the tornado 
became very strong and was stable for about 3.5 hours. 
The third flare seemed to have the biggest impact on the 
surrounding plasma or at least produced the most visible 
(at 195A) EUV wave. It also caused large active region 
loops to move towards the edge of the prominence cavity. 
Afterwards the loops swayed back and the cavity started 
to slowly (~ 2.5 km s _1 ) but visibly expand over the next 
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3.5 hours. During this time the prominence rose into the 
cavity creating the tornado. 

To understand the relationship between the flare, the 
prominence, the cavity and the tornado there are a num- 
ber of effects that need to be considered. The expansion of 
the cavity is clearly associated with the rise of the promi- 
nence plasma. One possible explanation for the expansion 
is that the loss of free magne tic energy from t he active 
region by flares and CMEs (Zhan g~fc Lowi[2QQ5l ) resulted 
in a contraction of the active region field, the Hudson ef- 



regK 

feet (|Hudsonll2QQQ[|Zhang fc LowBOOfllJanse fc Lowll2QQ/D . 
Subsequently the surrounding fields, including the cavity 
field, expanded to fill the vacated space. 

An idealized illustration of the Hudson effect, in terms of 
the force-free equilibrium between two bipolar fields before 
and after one of the fields has lost its free energy, is given 
in Fig. [51 Initially the active region field on the right has 
excess free energy, and the cavity field on the left is taken to 
be potential in order to simplify the calculation. During the 
eruption, the active region field loses free energy, decreasing 
its pressure since the magnetic energy density B 2 of a field 
B is also the magnetic pressure. This leads to an expansion 
of the cavity field to restore pressure balance in the system. 

We consider the 2D Cartesian domain < y < 7r, 
< z < oo, with x as an ignorable coordinate, taking 
the boundary to be a rigid perfect electrical conductor for 
simplicity. The initial state in Fig. [8^ is a continuous global 
solution for the magnetic field, B = ~Bi n iu a i of the force-free 
equations 



V x B = aB = 0; V • B = 0, 



(1) 



describing a bipolar (a = 0) potential field representing the 
cavity that occupies the partial domain < y < 27r/3,0 < 
z < oo. This field continues into a constant-a force- free field 
representing the active region that occupies the comple- 
mentary partial domain 27r/3 < y < 7r,0 < z < oo. Their 
interface y = 27r/3, is in force balance because the fields 
on its two sides exert equal magnetic pressure. Suppose the 
electric current of the constant-a force-free field is removed 
from the partial domain 27r/3<?/<7r,0<z<ooto rep- 
resent a flare-like loss of free energy. The distribution of the 
normal component of the field along the rigid boundary and 
the magnetic foot-points on the base z = cannot change. 
This boundary condition then determine the unique end- 
state, B en d, with a = everywhere in the domain (Fig.[Sb). 
The former bipolar field on the right has contracted down- 
ward and withdrawn into a partial domain of a finite height 
while the other bipolar field now occupies the whole infinite 
space above that height. The fluid-interface between the 
two fields is now an arc of finite length. These mathemati- 
cal solutions, given in the Appendix, serve to conceptually 
make specific a basic effect in the complex processes of a 
real energy release. 

In addition to the active region contraction, the 
CME/flare eruptions generated hydromagnetic distur- 
bances, seen as EUV waves, that perturbed the surround- 
ing low corona. They evidently caused swaying in the 
prominence and oscillations along the cavity boundary. 
These perturbations may have been intense enough to par- 
tially destabilize the cavity-prominence system, thereby 
contributing to or possibly even initializing the tornado- 
like activity. 

Our observational study supports the interpretation 
that the CME/flare activities in the active region had a 




Fig. 8. Idealized 2D 
Cartesian magnetic 
fields of the cavity 
and active region (a) 
before and (b) after 
the flare, in a domain 
< y < tt,0 < z < 
oo with z denoting 
coronal height. 



causal relationship with the nearby tornado. It is clear that 
the tornado was dynamically associated with the expan- 
sion of the prominence's helical field and its cavity. It is 
an open question as to whether the Hudson effect drove 
the cavity expansion and the helical field was responding 
dynamically to the expansion, or, the helical field loses its 
met a- stability as the result of hydromagnetic disturbances 
from the CME/flares from the active region. This physi- 
cal question is worthy of MHD time-dependent simulation. 
Evidence of the Hudson effect h as been reported in recent 
solar observational studies (e.g. ISun et al.l 120121 ; iLiu et aD 
2009). Our study shows the importance of this effect in 
relation to CMEs, flares and prominences, and provides 
observational motivation for theoretical MHD modeling to 
address the specific questions posed by this relationship. 
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Appendix A: Mathematical Illustration of the Hudson Effect 

We give the mathematical solutions for 2D magnetic fields in Fig. [8j which are straight forward to construct by standard 
techniques. Use the representation 



dA dA 

dz ' dy 



(A.l) 



in terms of the scalar flux- function A(y,z) and the component B x = Q(y,z) corresponding to the ignorable coordinate 
x. The initial field B iniUa i in Fig. 7a is given by A = A iniUa i and Q = Q initial where 

Ainitial = A pot = 2 sin exp(-j^); Q in itiai = 0, in < y < 2tt/3,0 < z < oo 

3 [27 

Ainitial = Afff = - sin (3y) exp(--z); Q initial = y y A in 2tt/3 < y < tt,0 < z < 00. (A. 2) 

The end-state field B en( ^ is the everywhere-potential field with A — A en d where 

00 

Aend = ^a n sinm/exp(-nz), (A. 3) 

n=l 

3 4 



1 2 x 3 4 27rn 1 
- and a n = sin — — 

3 7T 3 (4n 2 -9)(n 2 -9)' 



a 3 = - and a n = — - — sin^- — 2 - 2 — — , for n^3, (A.4) 



with Q = 0. Both ~B initial an d B en ^ have the same distribution of normal field component along the boundary of the 
domain < y < ir, < z < 00. 
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